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O ' Abstract 

<N 

We consider multi-dimensional assignment problems in a probabilistic setting. Our main 
results are: (i) A new efficient algorithm for the 3-dimensional planar problem, based on enu- 
' merating and selecting from a set of "alternating-path trees" ; (ii) A new efficient matching-based 

. algorithm for the 3-dimensional axial problem. 

1 Introduction 



CO 



o 



A (two-dimensional) assignment can be viewed as a set of pairs P = {(it,jt), t = 1,2, ... ,n} such 
that 

{il,h,...,i n } = {jl,32,---,jn} = N, (1) 



■ i.e., each row appears once in P, as does each column (and without loss of generality we may take 

(ii, . . . , i n ) = (1, . . . , n)). Given an n x n matrix of costs C = [Cij], the aim is to compute P that 
minimises C{P) = J2(ij)ep C«,j- This is a well-studied problem from the worst-case as well as the 
t^J- ■ probabilistic point of view and it is solvable in polynomial time. 

. In the standard probabilistic model for the assignment problem each entry Cjj is independently 

distributed as the exponential random variable with mean one, viz., Exp(l). There are numerous 
results related to this model: If Z n is the minimum value of C(P) then we have the remarkable 
result, conjectured by Parisi [13], 

n 

X'- E(Z n )=£±. 

This was proved by Linusson and Wastlund [11] and Nair, Prabhakar and Sharma [12]. Sec 
Wastlund [14] for a remarkably short proof of this. Earlier work proving that lirm^oo E(Z n ) = ("(2) 
was done by Aldous [1], [2]. 

In this paper we are concerned with the probabilistic analysis of multi-dimensional generali- 
sations of this problem. We consider two versions, Planar and Axial. Let us first consider the 
d-dimensional Planar model. Here we are given an n x • • • x n (i-dimensional matrix (tensor) C = 
[Cii,i 2 ,...,iJ> i- e -> a map C : [n] d — > R. An assignment is a set of n d-tuples T = {(i\,i2, ■ ■ ■ , i^), t = 
1,2, ... ,n} such that, in analogy with (1), for each dimension r G [d], {v>*ri • • • } = M- Geo- 
metrically, for r G [d], let an r-plane be a set of d-tuples of the form [n] r_1 x {x} x {n] d ~ r for some 
x € [n]. In the case of d = 2 a plane corresponds to a row or column of matrix C. An assignment 
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T is then a collection of d-tuples such that each of the dn planes contain exactly one d-tuple from 
T. 

The optimsation problem here is to compute an assignment T that minimises C(T) = 
z 2 i d )eT^h,i2,--;id- This P r oblem is NP-hard for d > 3 and the case of d = 3 is one of 
the original problems listed in Karp [10]. Some of its characteristics and applications are discussed 
in a recent book by Burkard, Dell'Amico and Martello [3]. Very little is known about the proba- 
bilistic behavior of the minimum of C(T) for d > 3. Grundel, Oliveira, Pasiliao and Pardalos 
[8] show that Z£ n — > whp in this case. At this point we can give some easy results on Z^ n which 
we state as 

Theorem 1 

(All proofs are given in the body of the paper.) The upper bound in this theorem is non-constructive, 
relying on recent work of Johansson, Kahn and Vu [9]. Our main result concerns a Bounded Depth 
Tree Search algorithm BDTS(fe). Here A; is a parameter that refers to the number of levels of 
search. It is unfortunate, but our approach only seems to give something interesting for d = 3. 

Theorem 2 Suppose that 1 < k < 7log 2 logn where 7 is any constant strictly less than 1/2. Then, 
whp 

(a) Algorithm BDTSffeJ runs in time 0(n 2k+2 ). 

(b) The cost of the set of triples T output by BDTS (k ) satisfies 

C(T) = 0{2 k n~ 1+9k logn) 

where 9 k = 2k +x^ ■ 

Note that for k such as |log 2 logn this is a "mildly exponential" running time, n°( logn \ yielding 
a solution which is an 0(n 1 / logn ) approximation to the optimum. After dealing with the Planar 
version, we will turn to the Axial version. Here we are again given an n x • • • x n d-dimensional matrix 
C. Geometrically, let a line be a set of d-tuples of the form x ■ ■ ■ {i r } x [n] x {v+2} x • • • {i^} 
for some r and ii, . . . , i r , i r +2 ■ ■ ■ ,id- In the case of d = 2 a line corresponds to a row or column of 
matrix C . An assignment T is then a collection of n d ~ l d-tuples such that each of the lines 
contains exactly one ci-tuple from T. 

The optimisation problem here is to compute an assignment T that minimises C(T) = 
S(n ii id)eT^h,i2„-,id- ^ ms P r oblem is NP-hard for d > 3 and this was proved in Frieze [6]. 
We will prove 

Theorem 3 The optimal solution value Z^ n satisfies the following: 

(a) Zf n = n(n d - 2 ) whp for d>3. 

(b) When d = 3 there is a polynomial time algorithm that finds a solution with cost Z where 

Z = O(nlogn) whp. 

This leaves the following open questions: 

PI What are the growth rates of EfZ^J and E[Z^J for d > 3? 
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P2 Are there asymptotically optimal, polynomial time algorithms for solving these problems when 
d > 3. 



P3 Frieze [5] gave a bilinear programming formulation of the 3-elimensional planar problem. There 
is a natural heuristic associated with this formulation (see appendix). What are its asymptotic 
properties? 

1.1 Structure of the paper 

We deal with the Planar version in Section 2. We start with the proof of Theorem 1 in Section 
2.1. Our next task is to analyse BDTS. We will analyse a three level version in Section 2.2. This 
provides intuition for the general case, analysed in Section 2.3, completing the proof of Theorem 2. 
The Axial problem is considered in Section 3. The lower bound in Theorem 3 is proved in Section 
3.1 and the upper bound in Section 3.2. 



Each term in the above sum is distributed as Exp(n d_1 ) and so has expectation l/n d_1 and variance 
l/n 2d ~ 2 . The Chebyshev inequality implies that the sum is concentrated around the mean. 

For the upper bound we use a recent result of Johansson, Kahn and Vu [9]. This implies that 
whp there is a solution that only uses d-tuples of weig ht at most ^SIL. The upper bound follows 
immediately. It should be noted that their proof is non-constructive. 

2.2 Two Level Version of BDTS 

In this section we consider a two level version of the algorithm BDTS. In this way we hope that 
to make it easier to understand the general version that is described in Section 2.3. With reference 
to Theorem 2, the two-level version means taking k = 3, 9 = 63 = 1/7. 
The heuristic has three phases: 

2.2.1 Greedy Phase 

The first phase is a simple greedy procedure. 
Greedy Phase 

1. Let m = n- n 1 - , J = K = [n], and T = 0. 1 

2. For i = 1, . . . , n\ do the following: 

• Let Cj iiifc = min{Cjj/ ;jt / : j' € J,k' € K}; 

• Add k) to T and remove j from J and k from K. 

1 We will often pretend that some expressions are integer. Formally, we should round up or down but it will not 
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matter. 
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At the end of this procedure the triples in T provide a partial assignment. Let 

(i,j,k)er 

Lemma 4 

2 

Z\ < —r-a wh P- 

Proof. We observe that if k) £ I then Ci j k is the minimum of (n — i + l) 2 independent copies 
of Exp(l) and is therefore distributed as Exp((n — i + l) 2 ). Furthermore, the random variables 
Cij,k, {hj,k) e T are independent. Using the facts that an Exp A random variable has mean 1/A 
and variance 1/ A 2 , 

"i, i rni+i ^ x j 

E(Zl) = ^ (n-z+1) 2 - J x=1 (n-x + l) 2 " n^' 



Now 
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Var(^) = £ J^TTW' ~ J™ = ° (E(Zl)2) 

1=1 

and the lemma follows from the Chebyshev inequality. □ 



2.2.2 Main Phase 

The aim of this phase is to increase the size of the partial assignment defined by T to n — 0(1). Let 
/ = I(T) be the set of first coordinates assigned in T, i.e., I = I(T) = {i : 3j, k s.t. k) € T}. 
Relabeling if necessary, without loss of generality we may assume that / = [\T\]. This phase will 
be split into rounds. We choose a small constant < a « 1 and let f3 = 1 — a. The aim of a 
round is to reduce the size of the set of unmatched first coordinates X(T) = [n] \ I(T) by a factor 
f3 while increasing the total cost of the matching only by an acceptably small amount. Thus we let 
x\ = n — n\ and xt = f3 l ~ l xi for t > 2. The aim of round t is to reduce |X(T)| from x% to xt+i- 
We continue this for to = l°gi//3( z i/ L) rounds where L is a large positive constant. Thus at the 
end of the Main Phase, if successful, we will have a partial assignment of size at least n — 2L. 

So suppose now that we are at the start of a round and that |X(T)| = x t . This is true for t = 1. 
Next let wq = 2n _12//7 logn and 

Wt = 2n- 6/7 x; 8/7 log 1/7 n for t > 1. 

At the start of each round we will refresh the array C with independent exponentials, at some cost. 
By this we mean that we replace C by a new array C where Cjj,fe < C[- k + wt-i and the entries 
of C are i.i.d. Exp(l)random variables. More precisely, suppose that during the previous round 
we determined the precise values for all Cjj,fe < wt-i and left our state of knowledge for the other 
Cjj,jfc as being at least wt-i- Then the memoryless property of exponentials means that 

,j,k - wt-i when C itj ^ > w t -i 

i j k — \ 

I fresh Ajj fc ~ Exp(l) otherwise 

has the claimed property. Thus we can start a round with a fresh matrix of independent exponentials 
at the expense of adding another wt-i to each cost. We note also that we can whp carry out the 
Greedy Phase only looking at those Cjj,jt of value less than wq. 
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Let T t denote the value of T at the start of round t and let I t = I(T t ),X t = X(I t ). In round 
t we will add At = [n — Xt + l,n — xt+i] to It- By relabeling if necessary we will assume that at 
the start of round t we have T = : 1 < i < n — x t }. To add i G A t to I t we find distinct 

indices j,k,p,q,r,s G It (distinctness is not strictly necessary) and replace 6 of the triples in I t by 
7 new triples: 

+ k) - - (k,k,k) + (j,p,q) + (k,r,s) - {p,p,p) - (q,q,q) - (r,r,r) - (s,s,s)+ 

(p, 6,6) + (q, 6, 6) + (r, 6, 6) + (s, 6, 6), (2) 

where 6, • • • ,6 are distinct members of X t , and each of the triples added in (2) is required to 
have (refreshed) cost at most wt- Roughly, we are assigning a new 1-coordinate i, this collides 
with previously used 2-coordinate j and 3-coordinate k, so the and (k,k,k) elements are 

removed from the existing assignment, 1-coordinates j and k are re-added as (j,p,q) and (k,r,s) 
thus colliding with the previous assignment elements (p,p,p), (q,q,q), (r,r,r), and (s,s,s), and 
finally 1-coordinates p, q, r, s are re-added as (p, 6,6) etc-, where the are elements not previously 
assigned. One may think of (2) as a binary tree version of an alternating-path construction; we 
will control the cost despite the tree's expansion. 

Putting Wt = wo + w\ + • • • + wt we see that if we can add one element to T at a cost of at most 
wt in refreshed costs, then in reality it costs us at most Wt; step (2) increases the cost by < 7Wt- 
Success in a round means doing this xt — xt+i times, in which case the additional cost of the Main 
Phase will be at most 7 times 

to to 

y^X x t - x t+iWt < x 1 (w + wi) + x t w t 

t=l t=2 

to 

< 3n~ 6 / 7 log n + 2x~ V V 6 / 7 log 1 / 7 n /T*/ 7 < An^' 7 log n. (3) 

t=2 

We must now show that whp it is possible to add xt — xt+i = ax t triples in round t with a 
(refreshed) cost of at most 7wt per triple. For this we fix t and drop the suffix t from all quantities 
that use it. We will treat refreshed costs as actual costs and drop the word "refreshed". 

We start by estimating the number of choices for assigning p. Ignoring other indices, the number 
of choices is distributed as the binomial Bm(u, 1 — e~ wx ) = Bin(z^, (1 — o(l))wx 2 ) where v = n — x. 

2 

Here 1 — e~ wx is the probability that for a given p, there exist 6,6 such that C p ^^ 2 < w. Note 
that 

wx 2 = 2(x/n) 6 / 7 log 1 / 7 n = o(l) and that wnx 2 3> logn 

and so the Chernoff bounds imply that, qs, 2 we can choose a set P of size exactly wnx 2 /2 = o(n), 
such that for each p G P there is at least one choice 6,6 £ X such that the triple (p, 6,6) is 
good, i.e., Cp ^ £ 2 ^ w - Given this set of choices P we find that the number of choices for q £ P 

2 

is distributed as the binomial Bin(^ — \P\, 1 — e~ wx ) and we can once again qs choose a set Q, 
disjoint from P such that \Q\ = wnx 2 /2 and each q G Q is in some good triple (q, 6,6) where 
6, 6 £ X- Similarly, we can choose sets R, S of choices for r, s, of size wnx 2 /2, such that P, Q, R, S 
are pairwise disjoint. 

Observation 5 Each £ G X is in Bin(xu, 1 — e~ w ) good triples of the form (p G P, £',£") and so 
qs it is in at most 

4n 1/7 log 1/7 n 
2wnx = —tz 

2 A sequence of events £„, n > are said to occur quite surely, qs, if Pr(£„) = 1 — 0{rC K ) for any constant K > 0. 
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such triples. 



We now discuss our choices for j and k. For a fixed j there are w 2 n 2 x 4 /4 pairs inPxQ and each 
has a probability f — e~ w of forming a good triple (j,p, q). Let j be useful if there is such a pair 
and useless otherwise. Then 

f 1 A \ 3 2 4 

w n x w n x 

Pr(j is useless) < exp < > < 1 . 

I 4 J 5 

It follows that the number of useful j £ Y = PliQU RU S dominates Bin(n — o(n), w 3 n 2 x 4: /5) and 
so qs we can choose a set J of useful j £ Y of size 

w 3 n 3 x 4 4ra 3 / 7 x 4 / 7 log 3//7 n 

= 3 = ° {n) - 

We can by a similar argument choose a set K of useful of this size disjoint from J and Y. 

Observation 6 A fixed p is in at most B'm(wn 2 x 2 /2, 1 — e~ w ) good triples (j,p,q) where (j,q) € 
J x Q and so qs every p is in at most w 2 n 2 x 2 such triples. 

Suppose then that in the middle of a round we have added y < ax triples to T. The number of 
£ G X that can be used in a good triple (p, £, 77) will have been reduced by y. The number of 77 will 
have been reduced by the same amount. It follows from Observation 5 that the number of choices 
for p will have been reduced by at most 2ax x 2wnx. By Observation 6 this reduces the number 
of choices for j by at most lax x 2wnx x w 2 n 2 x 2 + 7ax <C | J\ = u> 3 n 3 x 4 /6. The additional term 
+7ax accounts for the choices we lost because they have previously been used in this round. So 
our next % will get a choice of at least Bin((u> 3 n 3 x 4 /7) 2 , 1 — e~ w ) choices for a good triple k). 
So the expected number of choices is at least w 7 n 6 x 8 /49 = (2 7 /49)logn and then the probability 
there is no choice is o(n~ 1 ). This is sufficient to ensure that whp there is always at least one choice 
for every i. 



2.2.3 Final Phase 

We now have to add only 0(1) indices to /. At this point there is a problem with the bottom-up 
approach of the previous phase if x < 8, clearest in the case x = 1, say the single element n, when 
each of £1, ... ,^8 would have to be n, leading to an illegal assignment. Thus instead we will work 
top down. The details of this will cause more conditioning of the matrix, and therefore we refresh 
C after each increase in /, at an extra cost of w = Kn^l 1 log 1 / 7 n. So, if successful, the cost of 
this round is 0(W to + w) = 0(n~ 6 / 7 log 1/7 n). 
Let us now replace the notation of (2) by 

+ (hj,k)-(jij,j3)-(ki,k 2 ,k)+(j 1 ,p,q) + (k 1 ,r,s)-(p 1 ,p,p 3 )-(q 1 ,q 2 ,q)-(r 1 ,r, r 3 )-(s 1 ,s 2 ,s)+ 

(pi,h,P3) + (<?i,<72, J3) + {n,s 2 ,h) + (si, k 2 ,r 3 ), (4) 

where any subtracted triple such as (ji,j,j 3 ) denotes a previous match (we are no longer assuming 
the convention that such a triple would be (j, j,j)), and where i 2 , i 3 are unused 2- and 3-coordinates 
respectively. 

Fix j (and thus its previously matched companion indices ji,j 3 ) and let Zj be the num- 
ber of choices for p, q (with their previously matched companion indices pi,p 3 , qi, q 2 ) such that 
C(ji,p, q), C(pi, i 2 ,p 3 ), C(qi, q 2 , j 3 ) < w. This has the distribution Bi(B 2 (n,w)B 3 (n,w),w) where 
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Bi,B2,E>3 denote independent binomials, with B2 counting the good choices for p, B3 those for 
q, and B\ those for j using these p and q possibilities. Using Chernoff bounds on the binomials 
B 2 ,B 3 we see that whp Zj dominates B(n 2 w 2 /2, w) which dominates Be(n 2 w 3 /3), the Bernoulli 
random variable that is 1 with probability n 2 w 3 /3 and otherwise. The same holds for index k 
and (4) has been constructed so that choices for j, k are independent. So, the number of choices 
for j, k dominates Bin(n 2 , w(n 2 w s /3) 2 ) which has expectation fi(logn) and so is non-zero whp. 
This completes the analysis of BDTS when there are two levels. 



2.3 General 3-Dimensional Version 

We follow the same three phase strategy, k is a positive integer, 2 < k < 7 log log n. 



2.3.1 Greedy Phase 

This is much as before. Proceed as in Section 2.2.1 but taking 6 = 9k (recall 0's definition from 
Theorem 2) and defining n\ accordingly. Lemma 4 continues to hold. 



2.3.2 Main Phase 

Let 



a 



and let /?, to and Xt, t = 1, . . . , to be defined as in Section 2.2.2. Let It, X t ,A t have the same meaning 
as well. Now let wq = 2n~ 2 ^ x ~ dk) logn and 

w t = 2x^ 9k n 6k ~ 1 \og dk n for t > 1 

and 

W t = w + wi H \-wt = 



log 9fc n 



n 1 ~ 9k I 

The aim of round t is once again to add xt — xt+i new indices to It using triples with (refreshed) 
cost at most wt- We will assume that at the start of round t we have T = {(i, : 1 < i < n — x t }. 
In analogy with (2), to add i G A t to I t we will add 2 k+1 — 1 triples to T and remove 2 k+1 — 2 
triples, in which case the additional cost of the Main Phase will be at most 2 k+1 — 1 times 

to to 

y^(xt - x t+ i)w t < xi(wo + + x t w t 

t=l t=2 

to 

< 3rA^ log n + 2x~ dk n 6k - 1 log 9 " n ^ P~ 6kt < An^' 1 log n. (5) 

t=2 

The notation used in (2) is obviously insufficient. We imagine a rooted tree T of triples. The root 
will be p = (io, jo, ko) where io is the index to be added to It- The root is at level zero. The triples 
at odd levels are to be deleted from T and the vertices at even levels are to be added to T. Every 
triple at an odd level 21 — 1 will therefore have the form (p,p,p) where p € It- This triple will have 
one child (p, a, b) which will replace the parent triple in 1-plane p. li I < k then a,b € It and if / = k 
then a,b £ Xt- A triple u = (p, a, b) at an even level will have two children. By construction, u will 
be the unique triple in 1-plane p, but now we will have two triples in 2-plane a and 3-plane b. Thus 
the children of u are (a, a, a) and (6, b, b). This defines a tree corresponding to adding 2 k+1 — 1 and 
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removing 2 k+1 — 2 triples from T. We ensure that if u = (p, a, b) is a triple at an even level, then 
p,a,b do not appear anywhere else in the tree, except at the child of u as previously described. We 
do this so that additions in one part of the tree do not clash with additions in another part and 
then the additions and deletions give rise to a partial assignment. We also insist that if u = (p, a, b) 
is a triple at an even level then Cp^j, < w. We call such a tree feasible. We considered each level 
of r to be ordered so it makes sense to talk of the rth vertex of level 21 where 1 < r < 2 l . 

We now have to show that whp there is always at least one such tree T for each i G A t . We 
take the same bottom-up approach that we did in Section 2.2. We fix t and drop the suffix t from 
all quantities that use it. We start by estimating the number of choices for a p that can be in 
a triple (p, x, y) at level 2k. Ignoring other indices, the number of choices is again distributed 
as the binomial Bin(z^, 1 — e~ wx ) = Bin(z^, (1 — o{l))wx 2 ) where v = n — x = n — o(n). Note 
that wx 2 = K(x/n) l ~° k \og 6k n = o(l) and that wnx 2 = Q{n° k ) » logn. (Here our notation 
f(n) » g(n) means that f(n)/g(n) — > oo with n). So the Chernoff bounds imply that qs we can 
choose a set P of size exactly wnx 2 /2 = o(n), such that for each p 6 P there is at least one choice 
£i,£2 such that the triple (p, ^1,^2) is good, i.e., C p ^^ 2 < w. We will in fact be able to choose 2 k 
disjoint sets P^k, 1 < I < 2 k since replacing v by v — 2 k wnx 2 /2 will not significantly change the 
above calculations. (Here 2 k wnx 2 = 0(n l ~ 6k+6 k log efc+7 n) = o(n)). 

Observation 7 Each £ € X is in Bin(xi^, 1 — e~ w ) good triples of the form (p £ P^kiS,, ■) an d so 

- ° gw J 3> log n). 

Let 

uq = wnx 2 1 2 and i>i + i = wnvf /2 for < / < k. (6) 
The solution to this recurrence is 

2 !+l_l 



, 7 ^ (^) x 2l+1 = (nlognf l+1 -^ 2k+1 - 2l+1 ^ k 



Observe that v\ increases with /. Note also that if I < k — 2 then 

2 k 



< w fe 2 _ 2 = 2 (£) fe log( 2fc ^n = (l), (7) 

2 2 2 / i \ 29 k 

wnx / n log n x 



wnui > wnuQ = = 2 ^ J » logn. (8) 

We now have the basis for an inductive claim that qs if / < k — 1 and u = (p, a, b) is a triple at an 
even level 2(k — I) then there are at least v\ choices for p such that there exists a triple u = (p, a, b) 
with C u < w and a feasible tree T u with u as root and depth 21 + 1. Our analysis above has 
proved the base case of / = 0. Imagine now that we are filling in the possibilities for the rth triple 
(p,a,b) at level k — I. We fill in these possibilities level by level starting at level 2k. Imagine also 
that we have identified vi-\ choices for each of a, b. This can be an inductive assumption, so for 
example a will have to be a possible selection for the first component of the (2r — l)st triple at level 
2{k-{l-l)). 

For a fixed p, conditional on our having selected exactly v\_\ choices A, B for a, b, let p be useful 
if there is a pair (a, b) G A x B with C rha f, < w and useless otherwise. Then, using (7), 

r 9 -1 ^ wu i-i 
Pr(p is useless) < exp {— wv{_ 1 j < 1 — . 



8 



It follows that the number of useful p that have not been previously selected dominates Bin(n — 
o(n), 2wuf_ 1 /3). Here o(n) = Yl s <i 2 k ~ s wv 2 bound the number of forbidden p's. It follows that qs 
we can choose a set of useful p's of size wv 2 _ 1 /2 = o{n). We can do this so that each node of V gets 
distinct choices. 

Observation 8 A fixed a is in at most Bin(nz^_i/2, 1 — e~ w ) good triples (p, a, b) feasible for level 
2(k — I) and so qs every a is in at most wnv[_\ such triples, see (8). 

This completes our induction. We now apply the above to show that round t succeeds whp. 

Suppose that in the middle of a round we have added y < ax triples to T. The number of £ £ X 
that can be used in a good triple (p, £, rj) at level 2k will have been reduced by y. Thus the number 
of choices for p in any triple in this level will have been reduced by at most 2 k x 2 x ax x 2wnx, 
see Observation 7. This reduces the number of choices for p in a triple at level 2(k — 1) by at most 
2 k+2 awnx 2 x wuvq = 2 k+s awnv 2 , see Observation 8. So let \x,\ denote the number of choices for p 
in triples p(, ., .) at level 2(k — I) that are forbidden by choices further down the tree. We have just 
argued that fii < 2 k+3 awni>o- In general we can use Observation 8 to conservatively argue that 

m < wnvi_i(in_i + 2 k+1 ax). 

It follows that for I > 2 we have 

ill m_i 2 k+2 ax M-i / x \ ( 2 _1 ) efc u;_! ,,, 9 / x \ 0k 
tn < 2 ^i_J_ _^ < 2 h±J_ + 2 k+2 a — < 2^- + 2 k+2 a 1 ' 



v x vi_ x vi_ x " vi_ x " \n\ogn) v x _ x \nlogn 

It follows that 



A^l < ^-2^1 + 2 2k+l a ( " < 2 2k+2^ 

Vk-i v x \nlognJ 



We see that at the root there will still be at least (1 — 2 2k+2 a)vk-i choices for jo, ko- So i$ will 
get a choice of at least Bin((l — 2 2fc+2 a) 2 z/|_ 1 , 1 — e~ w ) choices for a good triple {io,jo, ko). So the 
expected number of choices is at least 2wv 2 _ l /3, our choice of a implies this. Now wv^._ 1 = 21ogn 
and this is sufficient to ensure that whp there is always at least one choice for every io- 



2.3.3 Final Phase 

We can execute the Main Phase so long as x > 2 k . Now assume that 1 < x < 2 k . We now 
have to add only 0(1) indices to /. This time we refresh C an 0(2 k ) number of times at an 
log^ fc n 

extra cost of w = — . „ each time we add an index. So, if successful, the cost of this round is 

0{W to +w)=o(^). 

We first make an inductive assumption: We have a partial assignment / where |/| < n — 2. (The 
reader might think that we should assume |/| < n — 1, but here we use the induction hypothesis 
after one more index has temporarily been deleted from /, prior to a replacement). Assume that the 
matrix C is unconditioned and i I: Then we can in 0(n 2e ) time whp find a set P of size (with 
x = 1 in definition (6)) and a collection Q p ,p € P of sets of size vg_\ such that for each (p, q £ P p ) 
there is an assignment P' with (i,p,q) G P' and I(P') ^ I(P) and C(P') = C(P)+C(i,p,q)+0(w). 
This is true for £ = 1 since we can make the changes 

+(hp,q) - (pi,p,P3) - (qi,Q2,q) + (pi,*2,P3) + (qi,Q2,h) 
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where 12,13 are unused 2- and 3-coordinates respectively. The number of choices for p, q are inde- 
pendent TMn(n,w). 

For the inductive step, we first refresh the matrix C. Then for each p G [n] we let I' = I — {p} 
and apply the induction hypothesis to generate vf_ 2 choices of assignment that add back p\ to 
We find that whp at least wuf_ 2 /2 = of these have C(pi, ., .) < w . Let this set be P. Now 
refresh C again and apply the same argument for each p £ P to generate choices Q p for p. This 
completes the induction. 

Now let t = k and refresh C one more time. Let P, Q p ,p G P be the sets of size Vk-i promised 
by the above argument. We have Bin(zv| _ v w) choices of j, k which can be used to add i ^ I to / 
at a cost of 0(w). In expectation this is 21ogn and so we succeed whp. 

For the execution time of the algorithm we simply bound the number of possible trees T. This 
completes the proof of Theorem 2. 

3 Multi-Dimensional Axial Version 

We turn to the proof of Theorem 3. 

3.1 Lower bound 

It is clear that Z^ n > Z\ + Z2 + • • • + Z n d-2 where Zi is the minimum cost of the 2-dimensional 
assignment with cost matrix Aj^ = Ci 1: ___ t i d _ 2 j : k- We know that Zj > (1 — o(l))C(2) whp and the 
Zi's are independent. It follows that whp Z^ > (1 - o(l))n d ~ 2 C(2) > 3n d ~ 2 /2. 

3.2 Upper bound for d = 3 

For the upper bound we need a result of Dyer, Frieze and McDiarmid [4]. We will not state it in 
full generality, instead we will tailor its statement to precisely what is needed. Suppose that we 
have a linear program 

P : Minimize (Fx subject to Ax = b, x > 0. 

Here A is an m x n matrix and the cost vector c = (c±,C2, ■ ■ ■ ,c n ) is a sequence of independent 
copies of Exp(l). Let Zp denote the minimum of this linear program. Note that Zp is a random 
variable. Next let y be any feasible solution to P. 

Theorem 9 ([4]) 

E(-Zp) < m max yj. (9) 
j=l,2,...,n 

Furthermore, Zp is at most 1 + o(l) times the RHS of (9), whp. 

Now consider the following greedy-type algorithm. We find a minimum 2-dimensional assignment 
for 1-plane i = 1, we then find a minimum assignment for 1-plane i = 2, consistent with choice for 
1-plane i = 1, and so on: 
Greedy 

1. For i = 1, . . . , n do the following: 

• Let G = K n>n \ (Mi U M 2 U • • ■ M;_i); 

• If (i, k) G E(G) let A jjk = dj tk . 

• Let Mi be a minimum cost matching of G using edge weights A. 
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The output, Mi, M2, . . . , M n defines a set of triples T = k) : (j, k) € Mj}. We claim that if 

Zi = A(Mi) then 

E < z .) s ^Vr. do) 

For this we apply Theorem 9 to the linear program 

Minimise Aj^Xj^ subject to 

(j,k)eE(G) 

xj,k = i, j = 1, 2 



5 5 1 



n 

fc: (j,k)€E(G) 

^2 x jtk = l, fe = l,2,...,n 

j: (j,k)eE(G) 

x j)k >0, j,k = 1,2,..., n. 

We note that there are 2n constraints and that x^fc = l/(n — i + 1) is a feasible solution. With 
Theorem 9, this implies (10) and the upper bound in Theorem 3 for the case d = 3. 
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A Bilinear Programming Formulation 

Frieze [5] re-formulated the 3-dimensional planar problem as 



i,j,k=l 

where Pa is the bipartite matching polyhedron Y17=i x i,j = 1 = Sj=i x i,ji f° r an 1 ^ h3 ^ n - 

Now denote the objective above by C(y,z). The following heuristic was used successfully in a 
practical situation [7]: 

1. Choose yo,zo arbitrarily; Z$ = C(yo,zo); i = 0. 

2. Repeat until Zi + \ = Z{. 

• Let yi+i maximise C(y,Zi). 

• Let Zi + i maximise C(yi + i,z). 

• Z i+ i = C(y i+1 ,z i+1 ). 

• % = i + 1. 



n 



Minimise 
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